Non-Fermi Liquid Behavior Induced by Resonant Diquark-pair Scattering in Heated 
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We show how the quasiparticle picture of quarks changes near but above the critical temperature 
T c of the color-superconducting phase transition in the heated quark matter. We demonstrate that 
a non-Fermi liquid behavior of the matter develops drastically when the diquark coupling constant 
is increased owing to the coupling of the quark with the pairing soft mode: We clarify that the 
depression and eventually the appearance of a gap structure in the spectral function as well as 
the anomalous quark dispersion relation of the quark can be understood in terms of the resonant 
scattering between the incident quark and a particle near the Fermi surface to make the pairing soft 
mode. 
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I. INTRODUCTION 

The recent data in the RHIC experiment suggest 
that the matter created by the RHIC seems to be 
a strongly coupled system with possible quasi-bound 
hadrons containedP,0. Furthermore some lattice calcu- 
lations, though in the quenched approximation, are con- 
sistent with the possible view that heavy-quark bound 
systems such as J/ip survive the deconfincment transi- 
tion at finite temperature T0. These developments have 
a possibility to change drastically the simple picture of 
the QGP phase that the system is composed of almost 
free quasi-particles, although such nontrivial properties 
of the QGP phase that it may contain quasi-hadronic 
excitations had been suggested earlier^, Q • 

In the present Letter, we shall argue that the dense 
quark matter at relatively low temperature can have also 
non-trivial properties, i.e., a non-Fermi liquid ones, if 
the system is close to the critical temperature T c of the 
color-superconducting phase transition 0, Q on a rather 
generic ground. 

First we notice that the low energy effective models 
of QCD show that the diquark gap at zero temperature 
may become as large as A ~ lOOMeV0, at lower 
densities such as those corresponding to it, ~ 400MeV 
with fx being the baryon-number chemical potential of 
the quark. Accordingly, the ratio r^ = A//i, which is 
a measure of the ratio of the inter-particle distance to 
the pair coherence length, may become as large as 0.2 to 
0.3. This value is much larger than those in the metal 
superconductors, i.e., ~ 0.001, which account for the 
validity of the mean-field approximation a la BCS the- 
ory for the electric superconductors. Thus one sees that 
large fluctuations of the diquark-pair are expected, which 
may invalidate the mean-field approximation, in the color 
superconductors 0, . 

In fact, our previous work with the fixed pairing cou- 
pling showed 0, that there exists the precursory soft 
mode composed of the diquark-pair field even well above 



T c . One of the remarkable points there was that the cou- 
pling of the quark with the soft mode not only modifies 
the quark dispersion relation but also causes a depression 
of the quark spectral function around the Fermi energy; 
the depression was so large that a "pseudogap" is formed 
in the density of states (DOS) of the quark even with the 
small diquark coupling that was employed. 

The above finding is interesting because they suggest 
that the properties of quarks at T above but close to 
T c are altered from the typical Fermi-liquid to the non- 
Fermi liquid ones owing to the coupling to the precursory 
fluctuations of the diquark-pair field. The mechanism to 
cause such a non-Fermi liquid behavior, however, might 
not have been clarified enough in the previous work. The 
purpose of the present Letter is to elucidate the mecha- 
nism to realize such a non-Fermi liquid behavior and ar- 
gue that it is a generic feature of the heated quark matter 
in the vicinity of T c of the color superconductivity. 

For this purpose, we examine what happens in the 
quark properties in detail if the pairing fluctuations be- 
come stronger by increasing the diquark coupling con- 
stant Gc; we remark that the limiting case can often 
clarify physics. Incidentally it may be noticed that Gc 
indeed could have other values than that in the previ- 
ous work without any contradiction with other principles 
and phenomenology. We calculate the dispersion relation 
and the spectral function of the quark in the T-matrix 
approximation|0 0, Q . 

We shall demonstrate the following for the first time: 
(1) The quark dispersion relation is so largely modi- 
fied with the increasing Gc that it becomes seemingly 
multiple- valued around the Fermi surface near but above 
T c in the strong coupling regime. (2) The quark spec- 
tral function gets to have a developed gap-like structure 
rather than a depression near but above T c for the larger 
Gc- We shall clarify that the above non- Fermi liquid be- 
haviors can be nicely understood in terms of the resonant 
scattering\uk between the quark and a particle near 
the Fermi surface to make the pairing soft mode. 



II. BRIEF SUMMARY OF FORMULATION 

To describe a system at relatively low T and p, it is ap- 
propriate to adopt a low-energy effective theory of QCD 
|l7l Il8j . Here, we consider the two-flavor quark mat- 
ter in the chiral limit with the four-Fermi quark-quark 
interaction^^ 

C c = G c ^(^75T 2 A^ T )(^ T «7 5 T 2 A^), (1) 



where C — 17270 is the charge conjugation operator, and 
T2 and Xa denote the antisymmetric flavor SU(2)f and 
color SU (3) c matrices, respectively. The quark chemical 
potentials p of each flavor and color are taken to be the 
same We introduce the three momentum cutoff A = 
65OMeV0. As for the diquark coupling constant Gc, 
we vary it in the range 3.11GeV~ 2 < G c < 4.35GeV" 2 , 
which is the similar range used in the references |sl 1^. IrH 

m 

In order to see the quasiparticle picture of the quarks, 
we examine the spectral function p Q and the dispersion 
relation of the quark lo — The following formula- 

tion is essentially a recapitulation of that given in \V'\ . 

The spectral function p is given by the imaginary part 
of the retarded Green function G R (k,Lo) as 



po(k,w) 



4tt 



Tr[7°ImG' R (fc,w) 



(2) 



The Green function G R (k,Lo) is decomposed into the 
quark and anti-quark parts; 



G R (k,u) 



(G R (k, u)A-(k) + G R (k, Lo)A + (k)) 7 ° 



A_(fc) 7 ° 



A+(fc)7° 



R-(k, u))+ir) R+(k,uj) + ii]' 



(3) 



with the projection operators A T (fc) = (1 ± 7°7 • fe)/2. 
The dispersion relations of the quarks and anti-quarks 



: (fe) are defined by the equations 
ReR T (k,oj) = 0, 



(4) 



respectively. A remark is in order here: Because w^(fe) is 
merely the solution of the real part but not the whole part 
of the inverse of the Green function R T (k,LO T ), Wip(fc) 
may not correspond to the peak position of the spectral 
function and hence also may not represent physical exci- 
tations when the imaginary part of the Green function is 
large. 

Our point in the calculation of the quark Green 
function G (k, lo) lies in incorporating the diquark-pair 
fluctuations in th e q uark self-energy in the T-matrix 
approximation|12t Il3l Il4| as is diagrammatically shown 
in Fig. ^ all the possible anomalous behaviors of the re- 
sults will originate from the fact that the pair fluctuations 
actually form the soft mode of the color-superconducting 
phase transition. 
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FIG. 1: Feynman diagrams representing the quark Green 
function in the T-matrix approximation. The thin lines rep- 
resent the free propagator Qo, while the bold ones the full 
propagator Q. 



In the Matsubara formalism, the self-energy of quarks 
in the imaginary time S(fe,w n ) is given by 



£(p,w n ) 



= 8T 



E 



d 3 k ~ 

—-^E(p + k,cj n + uj' m )g (k,io' m ), (5) 



with u> n = (2n+l)irT being the Matsubara frequency for 
fermions and Go(k,u>n) — [(^h+aOt k-'y] being the 
free quark propagator. The T-matrix for the quark-quark 
scattering S(k,u n ) reads 



E(k,v n ) = —G 



1 



c 



I 



(6) 



G c Q{k,v n y 

with the lowest order polarization function Q(k. v„)\l^ 
and v n — 2mtT. The T-matrix in the real time 
, E R (k,ui) = S(fe, v n )\i l y n ^ u]+ i n has all the information 
about properties of the pair fluctuations. In particular, 
the dynamical structure factor 



S(k,u) = - 



Im ~ R (k,iu) 



2ttG 2 c (1 



(7) 



represents the excitation probability of the pair field for 
each lo and k at finite T. It has a prominent peak at 
the origin near T c according to the softening of the pair 
fluctuations. 

The analytic continuation of S(fc, uj n ) to the real axis 
from the upper-half complex-energy plane gives the self- 
energy in the real time T, R (k,uj) = Y 1 {k,Lu n )\i Uri ^ iA]+ i 71 . 
Using the projection operators A T (fc), the self-energy 
is decomposed into the quark and anti-quark parts 
E T (fc,w) - Tr[£ fl (fc,w)A T (fe) 7 ]/2, and R T (k 7 uj) in 
Eq. © is now found to be R T (k,uj) = lo + \i =p k — 
T, T (k,oj). Notice that po(k,ui) around the Fermi energy 
is given almost solely by the quark part G R , and hence 
the corresponding part of the self-energy S_(fe, uS) is re- 
sponsible for the quasiparticle properties of the quark 
near the Fermi surface. For later convenience, we rewrite 
the explicit form of ImS_(fc, lo) given in jl3j slightly us- 
ing S(k,uj); 



TmE-(k,io) 



nGl 



d 3 q 



-S{k + q,u + E q - n) { 1 



3 



(l_e-(-+^-^)/ T )tanh 



Eg - P 

IT 



-(1 + c -(-+^-M)/T)j + {Eq ^ 



(8) 



with E q — \q\ 



III. NON-FERMI LIQUID BEHAVIOR AND 
RESONANT SCATTERING 



In this section, we show the spectral function po(k, lo) 
and the dispersion relation of quarks lo — w_(fc) to see 
the quasiparticle picture of quarks. As was promised in 
the Introduction, we vary the diquark coupling constant 
G c for a fixed e = (T - T c )/T c = 0.01. 

The upper panels of Fig. [21 show the spectral function 
Po(k, lo) for p — 400MeV and e — 0.01 with several values 
of Gc, while the lower panels the dispersion relation lo = 
L0-(k) together with the contour map of the strength of 
po(k,u>) represented by the shadow. 

The figures in the far left show those with the same Gc 
(T c = 40.04 MeV) as those presented in Refs. jHIH for 
comparison: Although the detailed properties of them 
are given in the previous paper[l3(, the essential points 
may be summarized as follows; (1) There appears a de- 
pression of the quasiparticle peak in the spectral function 
po(k,uj) around the Fermi energy u> = 0. This means that 
the life-time of the quasiparticles near the Fermi surface 
becomes short owing to the pair fluctuations. (2) The 
quark dispersion cj_(fc) shows a rapid increase around 
the Fermi momentum. 

In the middle and the far right panels of Fig. [21 Gc 
is increased 1.2 and 1.4 times larger than that in the 
far left panels, respectively; correspondingly, T c — > 59.82 
and 80.11 MeV. One sees from the figures that as Gc 
is increased, the depression of po(k,Lo) around the Fermi 
energy becomes more significant, and the quasiparticle 
peaks tend to be so completely depressed that the profile 
of pa(k, lo) around the Fermi momentum may be charac- 
terized more properly with the term of a gap-like struc- 
ture rather than with a depression. In the lower panels, 
lj- (k) around the Fermi-momentum shows more rapid 
increase for larger Gc, and eventually becomes triple- 
valued. Here one should notice that LO-(k) does not nec- 
essarily represent the real excitation spectrum since it is 
only a zero of the real part of the inverse of the Green 
function, as is noticed before. Indeed, one sees only two 
quasiparticle peaks in po(k,u>) for k — kp, correspond- 
ing to the first and third solutions of w_(fe). The reason 
of the disappearance of the middle solution of u_ (fc) in 
po(k,u>) will be given later. 

The above features of po(k,u>) andw_(fc) can be under- 
stood by the behavior of the quark self-energy £_(fc,w). 
In Fig. [3| we show the imaginary and real parts of 
£_(fc,u;) in the upper and lower panels, respectively, at 
k = kp for several Gc corresponding to each panel of 
Fig. 13 For each Gc, one sees a peak of |ImE_(fe,o;)| 



and a rapid increase of Re£_ (k, lo) around the Fermi en- 
ergy u) = 0, which is more significant for larger Gc- In 
fact, the size of Gc is important for the enhancement of 
£_ (fc, lo) because of the factor G c in Eq. © J2- We no- 
tice that the behavior of ReS_(fc, lo) can be understood 
by the growth of |ImE_(fc,w)| and the Kramers-Kronig 
relation 

Re£_(fc,w) = --P J du)'ImE_(k,u>')/(ui-ui'). (9) 

The peculiar behavior of w_(fc) can be understood as 
follows. First, recall that u>_(fc) at k — kp is the solu- 
tion of Re R_(kp,Lu) = lo — Re T,-(kp,u>) = 0. Accord- 
ingly, the solutions LO-(kp) are given graphically by the 
crossing points of ReE_(fc_p,w) and lo denoted by the 
straight dash-dotted line in the lower panel of Fig. [3] 
One sees how there eventually appear the three solu- 
tions of ui^(kp) for large Gc, as mentioned before. A 
remark is in order here: The group velocity as defined 
by v g = du^(k)/dk is seemingly larger than the speed 
of light near LO-(k) = 0. However, it corresponds to the 
peak position of the imaginary part of the self-energy 
|ImE_(fc, which means that the excitations around 
the origin have a quite large damping rate, and hence 
there does not appear a peak in the spectral function. 
The solution near u> = thus does not represent a phys- 
ical excitation spectrum of the quasiparticle. 

Now let us discuss the mechanism to realize the gap- 
like structure in the quark spectrum in Fig. [21 Due to the 
softening of the diquark-pair fluctuations near T c 
a particle near the Fermi energy is scattered by the soft 
mode and creates a hole, while a hole can create a parti- 
cle by absorbing the soft mode as shown in Fig.^Ja). We 
remark that the hole must be also created near the Fermi 
surface because of the energy conservation^^. More in- 
tuitively, the incident particle and a particle near the 
Fermi surface, which may be within the Fermi sea or ther- 
mally excited one, make a resonant scattering to form the 
pairing soft mode and vice versa. This resonating pro- 
cesses which are only effective around the Fermi surface 
induce a virtual mixing between the particles and holes, 
i.e. a (virtual) Bogoliubov transformation! In other 
words, because the particle and hole energies to — k — p 
and lo = p — k cross at the Fermi energy as shown in 
Fig. [5] the mixing between the particles and holes leads 
to the level repulsion of the energy spectrum of the parti- 
cle and hole, which makes the gap-like structure as shown 
in the right panels of Fig. [21 This mechanism which in- 
duces the gap-like structure owing to the resonant scat- 
tering to make the softening pairing mode is also known 
in the condensed matter physics [16|. 

We shall next show how the behavior of the quark spec- 
trum can be nicely reproduced by quantifying the notion 
of the resonant scattering. As was seen above, the soften- 
ing of the pair fluctuations near T c is responsible for the 
resonant scattering to become effective. Since the soft- 
ening is described as the enhancement of S(k,co) around 
the origin in the k-Lo plane, let us approximate the dy- 
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FIG. 2: The spectral function po(k, uS) (upper panels) and the dispersion relation of quarks u = ui-(k) (lower panel) for 
fi = 400MeV and e = (T — T c )/T c — 0.01. In the far left, middle and far right panels, the diquark coupling constant is taken 
Gc = 3.11GeV" 2 , 3.73(= 3.11 x 1.2)GeV~ 2 and 4.35(= 3.11 x 1.4)GeV" 2 , respectively. 
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FIG. 4: (a) The self-energy of quarks above T c . The dashed 
line denotes the collective mode of the pair fluctuations. The 
virtual mixing between the particles and holes is induced by 
the collective mode, (b) The self-energy in the BCS theory. 
The gap A induces the mixing. 



FIG. 3: The imaginary and real parts of the self-energy of 
quarks S_(fc F ,w) at k F = fi = 400MeV for G c = 3.11,3.73 
and 4.35. As Gc is increased, the peak of ImE-(fcF,w) as 
well as the increase of ReE_(fcF,^) around the Fermi energy 
uj — becomes significant. The dash-dotted line in the bot- 
tom panel denotes cu — k — fi. For larger Gc, there exist 
three solutions of Re R-(k,ui) = 0, giving the multi- valued 
dispersion relation. 



origin; 



ttG c 



(10) 



namical structure factor with the delta function at the 



The meaning of the choice of the pre-factor containing 
a A in front of the delta functions will be clear later; A 
is to be related with the imaginary part of the T matrix 
S and the boson distribution function 0, 0] . Needless 
to say, the dynamical structure factor above T c in reality 
never vanishes in the whole uj-k plane nor takes the form 
of Eq. H1U|I even at T = T c ; the fluctuations with finite 
lu and k make a width of the quasiparticles. However, 
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FIG. 5: The quasiparticle energy of quarks. 



this simple treatment will be found to clearly illustrate 
that the peculiar behavior of the quark spectrum can be 
nicely accounted for by the resonant scattering. 
A simple calculation with Eq. © leads to 



ImS_(fe,w) = -ttA 2 (5(w + |fe| - (j,), 



Re£_(fc,w) = A 2 



1 



Ifcl 



(11) 
(12) 



Then, i?_(fe,w) = gives real number quasiparticle en- 
ergies 



w_(fc) - ±V(|fej -m) 2 + A 2 



(13) 



which have the same form as that in the usual BCS theory 
with the energy gap A; notice that the quark propaga- 
tor for the quarks which interact with the pairing mode 
with the dynamical structure factor Eq. I|10[l can be di- 
agrammatically interpreted as Fig. 0Jb), i.e., the quarks 
interacting with a static source of the pair field A. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this Letter, motivated with the findings in our previ- 
ous papers[l2l[l3|]j we have tried to elucidate generically 
that the precursory pairing soft mode modifies the quasi- 
particle picture of the quark so much that the heated 
quark matter would become a non-Fermi liquid in the 
vicinity of T c of the color superconductivity. For this pur- 
pose, we have examined the changes of the quark proper- 
ties in detail when the diquark coupling is increased. We 
have found that in the strong coupling regime the quark 
dispersion relation is modified so greatly that it becomes 
multiple- valued around the Fermi surface near but above 
T c . The stronger coupling also changes the quark spec- 
tral function as to have a gap-like structure rather than a 
depression. We have shown that these non-Fermi liquid 



behaviors can be understood in terms of the resonant 
scattering between the incident particle and a particle 
around the Fermi energy to make the pairing soft mode. 

Although the present work is based on a model calcu- 
lation, the mechanism proposed here for realizing a non- 
Fermi liquid behavior of fermion systems can be model- 
independent. This is because the essential ingredient 
here is the scattering of a particle near Fermi surface to 
make the soft mode which is inherent to any second-order 
or weak first-order phase transition. It would be thus 
intriguing to explore the quasiparticle picture near the 
critical temperature of other phase transitions of QCD 
matter, including the chiral phase transition at finite 
temperature^- 

We should notice that there is another mechanism 
which causes a non-Fermi liquid behavior of ungapped 
quark matter: A resummed perturbation theory shows 
[21j that the coupling of the quark with unscreened long- 
range gauge (gluonic) fields gives rise to vanishing of both 
the residues of the quasiparticle excitations and the group 
velocity v g — dE/dp at the Fermi surface at small tem- 
perature without recourse to any phase transition. No- 
tice that the vanishing group velocity implies the infinite 
density of states of quasiparticles, which is highly in con- 
trast with our result, i.e., the pseudogap formation as a 
precursory phenomenon of the CSC at finite T. Although 
the perturbation theory adopted in [2l| is only applicable 
for extremely high-density quark matter, it would be in- 
teresting to see how the non-Fermi liquid behaviors due 
to unscreened gluonic fields survive and compete with 
those owing to the soft mode of CSC at moderate den- 
sity and temperature. 

In this work, we have employed the non-selfconsistent 
T-matrix approximation which is essentially a linear ap- 
proximation for the fluctuations. Nevertheless the results 
obtained in this approximation can be close to reality so 
lon g as e = (T—T c ) /T c is not so small as noted in |12lll3| : 
see |l5j. In the present work, the strong fluctuations of 
the pair-field were induced with the varied Gq for fixed 
e = 0.01. Similar results may be obtained for smaller e 
but with a fixed Gc since the fluctuations diverge as £ 
goes to zero. 

The situation, however, might not become so simple 
when the system may enter the Ginzburg region where 
the nonlinear fluctuation effects play an essential role and 
somewhat moderate behavior may be realized for the 
quark properties. A renormalization-group treatment, 
for instance, would be necessary to incorporate the non- 
linear effects of the fluctuations, which is beyond the 
scope of the present work and left for future investiga- 
tions. 
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